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Abstract 

We consider the regime in which the bands of the torsional acoustic (TA) and the hydrogen-bond- 
stretch (HBS) modes of the DNA interpenetrate each other. Within the framework of a model that 
accommodates the structure of the double helix, we find the three-wave interaction between the 
TA- and the HBS-modes, and show that microwave radiation could bring about torsional vibrations 
that could serve as a pump mode for maintaining the HBS-one. Rayleigh's threshold condition for 
the parametric resonance provides an estimate for the power density of the mw-field necessary for 
generating the HBS-mode. 
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I. INTRODUCTION 



It is generally accepted that the conformational dynamics of the DNA relies essentially 
on elastic vibrations of the DNA molecule in the region of 10 9 ± 10 12 Hz, According 
to Kim and Prohofsky, 

0,0 

, the latter includes the two domains, which correspond with 
different degrees of freedom of the molecule: (1) acoustic modes, which do not involve the 
hydrogen bonds; (2) modes that stretch the hydrogen bonds between the base-pairs (the 
HBS modes). Local minimum of the frequency is characteristic for the HBS-modes, Q], jjjl; 
its position depending on the choice of the band. The vibrations of the DNA were observed in 
the low-frequency Raman scattering, [4], Q], and the Fourier-transform infra-red absorption 
experiments, [(J, to a large extent depending on the state and composition of solvent. But, 
the experimental data, is not conclusive as to the relative positions of the TA- and 
HBS-modes. If the acoustic torsional bands interpenetrate the HBS-modes, the torsional 
vibrations of the double helix could change periodically the elastic constants of inter-strand 
motions, and thus provide a supply of energy for the HBS-modes. Then a torsional acoustic 
band, which has the double frequency with respect to that of the HBS-mode, could provide 
a means for maintaining an HBS-mode through the parametric resonance. 



II. THE ELASTIC DYNAMICS OF TORSIONAL AND INTER-STRAND MODES 

While considering the dynamics of the DNA one has to take into account: (1) the DNA 
having the two strands; (2) the base-pairs being linked by the hydrogen bonds; (3) the helical 
symmetry. We shall utilize a one-dimensional lattice model for the elastic properties of the 



DNA which accommodates t 
El Hasan and Calladine 



rese requirements. 
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set up the scheme for the internal geometry of the double 
helix of the DNA, which describes the relative position of one base with respect to the other 
in a Watson-Crick base-pair and also the positions of the two base-pairs. This is achieved by 
introducing local frames, or triads, for the bases and the base-pairs, and translation-slides 
along their long axes. We follow the guidelines of paper [3], but aiming at a qualitative 
description of the DNA dynamics use a simplified set of variables. We shall describe the 
relative position of the bases of a base-pair by means of the vector Y directed along the long 



axis ( y— axis of paper Q, see also paper [8]); Y being equal to zero when the base-pair is 
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at equilibrium. Thus, Y describes the base-base stretch corresponding to the displacement 
of the two-bases that make up the base-pair away from each other along the long axis y, 

The relative position of the base-pairs are described by the torsional angles 4> n , which give 
deviations from the standard equilibrium twist of the double helix. Thus a twist of the DNA 
molecule, which does not involve inter-strand motion or mutual displacements of the bases 
inside the pairs, is determined by the torsional angles <ft n that are the angles of rotation of 
the base-pairs about the axis of the double-helix. The twist energy of the molecule is given 
by the equation 



I • T 



in which / is the moment of inertia, and r and k are the twist coefficients, which for the sake 
of simplicity and taking into account the qualitative picture at which we aim, are assumed 
the same for all the base-pairs. 

Interstrand motions should correspond to the relative motion of the bases inside the 
base-pairs, therefore the kinetic energy due to this degree of freedom may be cast in the 
form 

n z 

where M is the effective mass of a couple. For each base-pair we have the reference frame in 
which (1) z-axis corresponds to the axis of the double helix, (2) y-axis to the long axis of the 
base-pair, (3) x-axis perpendicular to z- and y- axes (see Fig. 1 of paper [7]). At equilibrium 
the change in position of adjacent base-pairs is determined only by the twist angle Q of the 
double helix. We shall assume Q = 27r/10 as for the B-form of DNA. To determine the 
energy due to the inter-strand displacements we need to find the strain taking into account 
the constraint imposed by the helical structure of our system. For this end one may utilize 
the covariant derivative, as is done in paper jljj, but a simpler approach is possible. 

Let us confine ourself only to the torsional degrees of freedom of the double lattice and 
assume the vectors Y n being parallel to x-y plane, or two-dimensional. Consider the dis- 
placements Y n , Y n+ x for the two consecutive base-pairs, n, n+1. Since we must compare 
the two vectors in the same frame, we shall rotate the vector Y n+ i to the frame of the n-th 
base pair, 

r n+ f = R-\<p)Y n+1 

Here is the inverse matrix of the rotation of the n-th frame to the (n-fl)-one given 
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R{cf>) 



(1) 



by the equation 

cos <p — sm ■ 
sin cos ( 

The matrix R is 2 by 2 since the vectors Y n are effectively two-dimensional. Then the strain 
caused by the displacements of the base-pairs is determined by the difference 

\F back xr 
I n+1 ~ 1 n 

For this argument I am indebted to D.I.Tchertov. 

It is important that the angle <fi is given by the twist angle, Q, describing the double 
helix, in conjunction with the torsional angles ra , so that 

= + (f) n+1 - (f) n 

Therefore, the energy Tiy due to the base-base stretch, or the inter-strand stress, reads 



2 n 



'H) y A — y„ +— n- ] io--<p n+1 -<p n )Y n+l -Y n 

in which K and a are the torsional elastic constant and the inter-pairs distance, corre- 
spondingly. In summations given above n is the number of a site corresponding to the n-th 
base-pair, and n = 1,2,. . . ,N, N being the number of pairs in the segment of the DNA 
under consideration. The constant K is assumed to be the same for all the base-pairs. In 
fact, it differs from the situation in real life, but we are aiming at a rough qualitative picture. 
Besides it is possible to manufacture artificial double helices of the form (polyX).(polyY), 
in which the elastic constants shall not depend on the choice of the base-pair. The last 
term, e/2Y 2 accommodates the energy of the inter-strand separation due to the slides of the 
bases inside the base-pairs. The form of TCy corresponds with the fact that the equilibrium 
position of the double helix is the twisted one determined by Q, all <j) n being equal to zero. 
We suppose that the size of DNA molecule is small enough that it can be visualized as a 
straight double helix, that is not larger than the persistence length. Hence the number of 
base-pairs, iV < 150, approximately. The total energy TC e i as tic of the DNA molecule reads 

^elastic = T~i<t> + 7~(-Y (2) 

Preserving only terms up to the third order in (f> n , and Y n , we may transform TCeiastic 
given by Eq.(J2J), so that it takes on the form 



Ti 



elastic 



E 



1 12 r 



( M 4 2 r 



^ 2 n 



K , 



R-\Q)Y n+1 xY n 



(3) 



We have used the fact that the axis of the double-helix is directed along Oz-axis. 
Let us simplify the equation for TL e iastic with the help of the unitary transformation 

Y n = -j= On + iu l) , Y n = ^= On + <4) 



Note that the vectors Y n and u n are effectively two-dimensional, their third coordinates 
being equal to zero. On applying the Fourier transform 

1 



fn 



E 



-vnaq 



where q — j^m, m — 0, ±1, . . . , ±y we cast the equation for the energy Ti. e i as tic in the 



form 



n 



elastic 



= E 
+ E 



\k <fig 
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2 Q — aq 1*1 - £1 + aq 2 * 2 
sm — - — u q u q + sm — - — u q u q 



-iaq 



N 



-e iQ u\, u ql+q „ +e in u 2 q , u„, 



q'+q" 



(4) 



It is important that after the Fourier transform the variables u n verify the following equations 
for their complex conjugates 



*i 



iu 



IU 



~q "">-qi ">q "^-q (5) 

The interaction term in the last equation corresponds to the three-wave process, and may 
result in resonance. We shall utilize the fact for deriving the parametric maintenance of the 
u q modes, i.e. the HBS modes, (see below). 

One can obtain, in the usual way, the equations of motion for u q , a = 1, 2 and <p q , from 
the equation for the energy indicated above. The essential point is the effects of dissipation, 
which are due to ions in the close neighborhood of the molecule and water effects, see jl^ . 
The dissipation can be accommodated by writing down terms linear in u q , <p q . We shall take 
into account external force, or torque T q , only in the equation for <p q , for it corresponds to 
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external degrees of freedom of our model. Thus, the equations of motion can be cast in the 
form 



4i^sinf2 ^ • » „ , . 

+ T= Y e" q </y< a > = 0, 6 

Ma 2 VN^ Vq q ' q y 1 

+ ' IcfiJN \ Uq ' Uq ' ql = q (7) 



in which the dispersion laws for the fields u q , a = 1, 2, and 4> q are given by the equations 

^ 2 - AK sin 2 M + t- 1 ) ^ + 6 ^ 

^ " Ma 2 " Sm 2 + M' (8) 

W 4 = 777 sm ^7 ( 9 ) 



We see that the spectrum of <fi q has a typical acoustic character, whereas that for u q has 
a local minimum determined by the helical twist, Q. Thus, the spectrum of our model 
is in qualitative agreement with conclusions of 

a, a. 

The nature of the torque is to be 
specified elsewhere (see Section IIIIJ) . For the moment, we shall consider general dynamical 
phenomena to which the torque may be conducive. 

Let us suppose that for one thing the amplitudes of the HBS-modes given by u q be 
so small that the quadratic term in Eq.(J2J) can be neglected, and for another the external 
torque T q be appreciable enough to maintain the vibration of the torsional mode <p q . Thus, 
we can visualize the latter as a pump mode that interacts with the HBS-mode u q through 
the non-linearity in Eq. ([6|) . Since we aim at studying the problem of mw-radiation (see 
next section), we shall confine ourself to the case of the torque T q being non-zero only at 
q = q* = Q/a, and having the frequency 2u. Therefore, the forced wave, or the pump wave 
for the HBS-mode, has the form 

<P q , = e i2 ^ $ V , = e' i2ut $* &_ M . (10) 

To obtain larger values for the pump wave, (fi q , the resonance condition 

should be verified. The equations of motion for in the pumping regime read 
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« + + + 777777^7 ( 



Ma 2 ^/iV 



here 
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Note that the momentum conservation in the q-values is preserved, as required by the three- 
wave interaction. The equation indicated above can be cast in the matrix form 



1C + )u. 



(11) 



ii a + u 2 a u a + 7u il Q = (e l2ut K + e~ l2wt 
where JC and /C + are hermitian conjugate, and 

KL + fC — X ( — — ' ; — ) I A I 2 , Xti = 5a 
\Ma 2 v/jV / 

It is worth noting that Eq. ffTTj) is a kind of matrix Mathieu equation. In fact, we can apply 
to it Rayleigh's method for studying parametric resonance, [lOj. For this end let us look for 
the solution to Eq. (jll)) in the form of a series 



u{t) = A x e lwt + B x e' luJt + A 3 e liuJt + B 3 e~ l6u>t + ... 

On substituting the expression given above into Eq.([ll|) and preserving only the terms 
corresponding to e ±luJt , we obtain the equations 

-uo 2 + i lu uo) l + co 2 a ] A x + JCBx=0 
-uo 2 - i lu uo) X + uo 2 a ] Bx + /C + Ax = 

The compatibility condition of the equations indicated above can be cast in the form of 
determinant for the block matrix 



det 



uj 2 a -u 2 + Z7„ LO 



ujI-uj 2 - i>y u uo 







(12) 



Here Cj 2 is the matrix of frequencies given by Eq.fjHJ), and uo 2 and j u u are the scalar ones. 
Using the fact that for the range of frequencies we are considering, the matrix 

Qj 2 a - J 2 + uj 



is not degenerate, we may cast Eq. (jl2j) in a more amenable form given by the equation 



2K sin ft N 2 



(u4 -lo 2 -i 1u lo){lo 2 _ g -u 2 + i lu u)- [— ^= |A| 2 = (13) 



which is quite similar to the usual condition for parametric resonance. Solutions to Eq. (jl3|) 
are generally complex and therefore correspond to attenuated regimes. But there is a specific 
wave number, q res , for which the solution gives the real frequency u, and it is easy to see 
that it should satisfy the constraint 

Ula-q* = ^Iqi Q = Ires (14) 



Thus, we may cast the condition for parametric resonance in the familiar form, 



110]. 



(co 2 -co 2 ) 2 + 7 2 ^ 2 - f^^£V Ml" U r )) 
aq res J i \Ma 2 v^V / 



III. MICRO- WAVE IRRADIATION AND THE HBS-MODES 

We may use the results of the previous section for assessing the action of mw-radiation 
on the molecule of the DNA. The key point is accommodating the fact that the wavelength 
of radiation is by many orders of magnitudes larger than the characteristic size of the region 
of the molecule involved in the process. It was Chun-Ting Zhang, Q], who suggested a 
mechanism to overcome this difficulty The main point of Zhang's argument is that the 
helical configuration of the electric dipoles corresponding with the base-pairs makes the 
interaction U = —P ■ E of the dipole P and the field E dependent on angle, and therefore, 
different torsional momenta are applied at the base-pairs. For our model the equation for 
the energy of interaction between the dipoles of DNA and an incident micro-wave reads 

- E ■ R(nn + n )P o 

n 

Here R(nfl + <f> n ) is the rotation matrix given by Eq.flTJ), and P D is the dipole at site n = 0. 
Consequently, even though at the molecular scale the radiation has a plane wave configura- 
tion, it still twists the DNA molecule about the axis of the double-helix. Since the momenta 
changes periodically in time with the incident wave, the irradiation results in a periodic 
stress that may produce elastic vibrations in the DNA molecule. Zhang suggested that the 
force may generate resonance vibrations, resulting in a cross-over mechanism which takes 
up initial torsion excitations and transforms them into longitudinal acoustic vibrations. 
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In the present paper we will try to combine Zhang's mechanism, ^J, and the excitations 
of the double-helix studied by Prohofsky and Kim, [4], 0, with the view of generating inter- 
strand waves in the DNA by mw-irradiation. In contrast to the original idea by Zhang, we 
do not utilize a cross-over into longitudinal acoustic vibrations, but employ the interaction 
between torsional oscillations and the inter-strand ones, i.e. the three- wave, given by Eq.(J3J. 

The main point is that by expanding the rotation matrix R(nQ + <f> n ) in the angles <p n 
and keeping only the first order terms, we may cast Zhang's interaction in the form 

Hz = ~ E K (E x P n ) 3 + const, P n = R{nty P Q (16) 

n 

in which P Q is the dipole vector at site n = 0. Next, by using Eq.(JJ) for the matrix R(nQ) 
and neglecting the constant term we may cast Eq. (jl6|) in the form 

W z = \ £ <t> n {e mn \{E x P ) 3 - i(E ■ P )} + e~ tnn [(E x P ) 3 + i(E ■ P )\ } 
On applying the Fourier transform for the (f) n , and utilizing the equation 

n 

we obtain the following expression for Zhang's interaction 



^z = y Y~ {0 9= £ {(E x P ) 3 - i(E ■ P )) + cj> q= _n a [(E x P ) 3 + i{E ■ P )]} 
Hence, the torque T q in Eq.© corresponding to Tiz is given by the equation 

Z Z* 

T = — 6 q _ g . + — 8 m * q if = VL/a (17) 

in which 

Z = ^- [{E*P )z + i{E-P ) 

It should be noted that ±g* are the local minima of the HBS-modes; q res given by Eq.(jl 
reads 

3 

q re s = -q* (18) 

It is worth noting that the wave numbers q* and q res correspond to the wavelengths of one 
and I turns of the double-helix. 
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Equations given above provide an opportunity for making numerical, order of magnitude, 
estimates, which enable us to assess the effect of mw-radiation on the HBS-modes. From 
Eq.(fT7j) we infer that the torque T has the size 



T oc e 2iut E P 

where E and P are the external field and the dipole moment of the base-pair, respectfully. 
Next, suppose that the resonance condition 

Q 

UJ q = 2lu, Q = Q* = — 

a 

be true, so that there is the resonance of the incident mw-radiation and the acoustic torsional 
mode <p qt . Therefore, we may expect that the action of the radiation on the torsional modes 
should be the largest possible. Then the amplitude of the pumping wave, (p qt , according to 
Eq.©, is of the order 



^ VN EP 

$ oc 19 

Next, we turn to Rayleigh's condition for the parametric resonance of the HBS-mode given 
by Eq. (fTK|) . For the pumping wave corresponding to Eq. (fTT)|) . it gives 

2 

Ma 2 Pfyu ) 
Hence we have the threshold condition 



/ 9 99 , f K sintt EP \ ' 

- <J +7^ 2 ~4 



2K sinVl EP , . 

7 ^-M^— (20) 

We suppose that the frequency of the HBS-modes be determined by the first term in Eq.fjHJ), 
so that the first factor in Eq. (j2(J|) does not differ much from unity. It signifies that the 
energies of the inter-strand separation per base-pair is approximately of the same order as 
that due to the twist of the relative positions of the two adjacent base-pairs. If so, we could 
have the estimate for the dissipative constants, at least by orders of magnitude, 

iui4><-r ( 21 ) 



On utilizing the expression 



S = -j-E x H, 
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for Pointing's vector, we cast the estimate given by Eq.(|21|) in the form 

7u7,<2j ] j^ (22) 

in which S is the power density of the interaction and c the light velocity If we assume 
oc 1 Debye and the inertia coefficient / oc 10~ 36 gr cm 2 , corresponding to the mass of the 
base-pair oc 10~ 22 gr, and the size oc 10 A, then for the power density S oc 100 mW/cm 2 , we 
have 

7«7<^<10 16 ^ 2 

The estimate suggests that the effect produced by mw-radiation is to be looked for at the edge 
of the GHz zone, for in this case the requirement on the line- width is less stringent. It should 
be noted that the crucial point in assessing the feasibility of experiments on mw-irradiation 
of the DNA, and its possible influence, is the part played by ions in ambient solvent. In 
fact, the irradiation may result in just heating the solvent, so that the dissipation due to 
the ions takes up all effects on the molecules of DNA. Generally, the thin boundary layer of 
water and ions close to the DNA-molecule may have an important bearing on the dynamics 
initiated by the incident mw-radiation and result in the overdamping of the molecule's 
torsional oscillations. Davis and VanZhandt, 12], put forward arguments that nonetheless 
the influence of the boundary layer could be effectively small so as to allow the survival of 
the effect caused by the mw-irradiation. But, so far there have been no definite arguments 
in this respect, and studying the mw-effect may turn out to be helpful for its understanding. 



IV. CONCLUSIONS : IMPLICATIONS FOR THE DNA BIOPHYSICS 

The elastic dynamics of the double helix could have enough structure for providing a 
means for stretching the hydrogen bonds of the base-pairs of DNA, or generating the HBS- 
modes. If the vibrational modes of the DNA are not overdamped by the ambient solvent, and 
the balance between energies supplied and dissipated is favorable, the "optimistic" forecast 
is that sufficiently developed resonance instabilities could bring about the breaking of the 
relatively weak H-bonds ( with free energy of less than 2ksT , see ^|), and even result in 
the denaturation of the molecule. On the contrary, if the dissipation is strong, which is quite 
likely, the maintenance of the HBS-modes could be expected only at the edge of the HBS- 
zone and for sufficiently strong pumping. It is important that in the range of frequencies less 
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than 10GHz, the absorption due to water is small, whereas close to 100GHz the permitivity 
of water changes drastically. 

The best technique for studying the H-bond stretching still remains the Raman spec- 
troscopy on which certain improvements have been made (see 3] and references therein). 
The choice of specific means for generating torsional excitations of the DNA is important 
and interesting. In this paper we have envisaged mw-irradiation of the DNA. In case the 
interpenetration of the acoustic and the HBS-modes takes place, we come to the conclusion 
that mw-radiation could maintain the HBS-modes, if the power density is sufficiently large, 
100 mW/ cm 2 or more. It is important that there is no need for long exposures of the sample 
to the radiation. Thus, the double helix could be a kind of mw-parametric amplifier for 
the HBS-modes, which are always present due to heat fluctuations; the torsional vibrations 
serving a pump beam for the HBS-modes. It is worthwhile to take into account that the 
parametric resonance can take place also for the frequencies p of the pumping wave that 
verify the equation p = 2u/n, n = 1,2,3.... Therefore, the constraint imposed on the 
frequencies of the AT- and HBS-modes is less stringent than it appears at first sight. If 
the effect be sufficiently pronounced, it may result in the formation of the bubbles of bro- 
ken H-bonds. In this respect it is worth noting that our estimate for the critical power 
density, 100 mW/ cm 2 , is by orders of magnitude larger than that officially prescribed, i.e. 
0.2 - lmW/cm 2 . 

The whole range of problems related to the direct action of mw-radiation on the DNA 
is highly controversial (see H. Lai's discussion of the subject, [16]). For one thing there 
is enough evidence for the conclusion that it has a bearing upon the functioning of living 
organisms; for another, "micro- wave" effect, considered as the resonance absorption of ra- 
diation, is difficult to distinguish from accompanying phenomena, for example heating, as 
well as to assess its real implications. In this respect utilizing the HBS-mode, which can 
be maintained by short pulses of radiation and detected and registered by the experimental 
methods currently employed in condensed matter physics, could be a valuable probe into 
the action of mw-irradiation on the DNA molecule. 
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